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Abstract. An abstract should be given 

e display by numerical calculation how rest frame spectral lines appear in the observed 
spectrum of gamma-ray bursts due to the Doppler effect in the fireball framework. The 
analysis shows that: a) in the spectrum of a relativistically expanding fireball, all rest 
frame lines would shift to higher energy bands and would be significantly smoothed; b) 
rest frame weak narrow emission lines as well as narrow absorption lines and absorption 
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line forests would be smoothed and would hardly be detectable; c) the features of rest 

^ , frame broad emission lines as well as both strong and weak broad absorption lines would 
•i-H . 

■ remain almost unchanged and therefore would be easier to detect; d) deep gaps caused 
$_i ' 

C$ ' by rest frame broad absorption lines would be significantly filled; e) a rest frame emission 
line forest would form a single broad line feature; f) the observed relative width of the 
rest frame very narrow line would approach 0.162; g) when the Lorentz factor T is large 
enough, the observed line frequency vu ne and the rest frame line frequency vo,ime would 
be related by vu ne w 2rVo,hne- We also investigate the effect of time dependence of the 
line intensity and the effect of variation of T. We find that the feature of rest frame dim- 
ming narrow emission lines would disappear when F is very large. The form of emission 
lines would be sharp on both edges when T varies with time. This phenomenon depends 
not only on the initial Forentz factor but also on the observation time. 

Key words, gamma-rays: bursts — gamma-rays: theory — radiation mechanisms: non- 
thermal — relativity 
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1. Introduction 



Due to the observed great output rate of radiation of gamma-ray bursts (GRBs), most 
models of the objects envision an expanding fireball (see e.g., Goodman 1986; Paczynski 
1986). The gamma-ray emission would arise after the fireball becomes optically thin, in 
shocks produced when the ejecta collide with an external medium or in shocks occurred 
within a relativistic internal wind (Rees & Meszaros 1992, 1994; Meszaros & Rees 1993, 
1994; Katz 1994; Paczynski & Xu 1994; Sari et al. 1996). As the expanding motion of the 
outer shell of the fireball would be relativistic, the Doppler effect must be at work. As 
pointed out by Krolik & Pier (1991), relativistic bulk motion of the gamma-ray-emitting 
plasma can account for some phenomena of bursts. However, in many cases, the whole 
fireball surface should be considered, and accounting for the effect the fireball surface 
itself would play a role (Meszaros & Rees 1998; Qin 2002, hereafter Paper I). 

As revealed by the observation of absorption lines of afterglows (Metzger et al. 1997), 
GRBs are confirmed to be events occurring within the environment of stars. Based on 
the assumption that there might be a strong magnitude field within the fireball and that 
the expansion would be relativistic, it was believed that synchrotron radiation would 
become a dominate mechanism (Liang et al. 1983). Unfortunately, the spectra observed 
so far are so different that none of the mechanisms proposed can account for most of the 
spectral data of the objects (Band et al. 1993; Schaefer et al. 1994; Preece et al. 1998, 
2000). 

Based on the assumption of the environment of stars, mechanisms of formation of 
emission as well as absorption lines for GRBs have been well established. Meszaros 
and Rees (1998) suggested that gamma-ray burst outflows may entrain small blobs or 
filaments of dense, highly ionized metal-rich material, which might be accelerated by 
the flow to Lorentz factors in the range 10 — 100. In the event of neutron-star colli- 
sions or black hole-neutron star collisions, neutrino-antineutrino annihilation can pro- 
duce electron-positron pairs. Thus it would be natural if there exist in the outer shell 
of fireballs some high energy emission lines such as the 6AkeV line and the bllkeV 
annihilation line (Ramaty & Meszaros 1981). As reviewed by Piran (1999), both absorp- 
tion and emission features have been reported by various experiments prior to BATSE. 
Absorption lines in the 20 — 40keV range were observed by several experiments. GINGA 
discovered several cases of lines with harmonic structure (Murakami et al. 1988; Fenimore 
et al. 1988). These lines were interpreted as cyclotron lines (reflecting a magnetic field of 
« 10 12 Gauss). Emission features near AOOkeV were claimed in other bursts (Mazets et 
al. 1980). However, so far, BATSE has not found any of the spectral features (absorption 
or emission lines) reported by earlier satellites (Palmer ct al. 1994; Band et al. 1996). 
Meszaros and Rees (1998) pointed out that within the relativistic fireball model an ob- 
served broadened spectral line might be a blue-shifted iron X-ray line. We suspect that, 
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in detecting line features, the Dopplcr effect associated with the expanding fireball might 
play a role. 

In the following we will investigate what one should expect in the spectrum if there 
are emission as well as absorption lines in the rest frame of fireballs. 



2. Effect on narrow emission lines 

We consider a fireball expanding with a definite value of Lorentz factor T and study only 
the core content of the Doppler effect in the fireball framework with the cosmological as 
well as other effects being temporarily ignored. The rest frame radiation is assumed to 
be constant and independent of direction. 

Let the fireball be observed at time t and at frequency v. Then the expected spectrum 
would be (Paper I) 



2vR 2 {t) V r' 2 h. v {t ,e,vo. e ) cos 9 sintf 

Jo (1-/JCOS*)* (1) 



D 2 r 3 

with 



R(t) = f3[c(t - t c ) - D] + R c , (2) 

where D is the distance of the fireball to the observer; 9 is the angle to the line of sight; 
vo,e is the rest frame emission frequency of the differential surface, of the fireball, with 
9; to,e is the proper emission time of the differential surface; t c is any coordinate time 
concerned; R c is the radius of the fireball at time t c ; and Io,v(to,e, vq.b) is the rest frame 
intensity of the differential surface. Frequencies Vq^ and v are related by the Dopplcr 
effect. The proper time to,0 and observation time t can be well linked by considering the 
travelling of light from the fireball to the observer (see Paper I) . 

Assume that during some period the radiation of the fireball is dominated by a certain 
mechanism and within this interval of time the rest frame radiation intensity can be 
expressed as: 

Io,v(to,e, v ,e) = hihfi)go. v {vQfi), (3) 

where go lV (vo,e) describes the dominant radiation mechanism while Io(io,e) represents 
the development of the intensity magnitude. For constant radiations, Jo(io,e) = ^o- In 
this case, the spectrum would be 

, 2irI R 2 (t)v r /2 , . cos6»sin6i in 

here the emission is assumed to cover the whole energy range. 

As mentioned above, there is no single mechanism proposed that can well represent 
all the observed spectra of GRBs. In practice, an empirical form called the GRB model 
(Band et al. 1993) was frequently, and rather successfully, employed to fit most burst 
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Fig. 1. The expected spectrum of a fireball with its rest frame radiation being the GRB 
form of a = —1 and fto — —2.25 plus the two emission lines defined by (8) and (9), 
observed at time t, where we take 2-KloVQ^ p R 2 (i) /D 2 = 1. The solid lines from the bottom 
to the top correspond to T = 1, 2, 10, 100, 1000, 10000, respectively. 

spectra (Ford et al. 1995; Preece et al. 2000). The relative intensity of the GRB model 
is (Band et al. 1993) 

9o,„cM = { „ Pi (2 + a °' G) -J 2 +-, G ) (5) 

)a0iO _* iO exp(/3o g _ Qo g)( ^ ) i +/3o , g ( ai > s^dks), 

where ao,G an d po,G are the lower and higher indexes, respectively. As an illustration, 
we shall adopt in the following the GRB form with the typical values of the parameters: 
a G = -1 and /3 ,g = -2.25 (Preece et al. 1998, 2000). 

To show the effect on both emission and absorption lines, we consider the following 
rest frame radiation: 

n 

9oA v o,e) = max{5o,„,G(fo,8) + e ,^(^o,e), 0} (6) 
where eo,u,i{vo,e) describes the relative intensity of lines: 

^ M = exp(2 + ao ) CXp [ 23710^ ] (?) 

Now we consider two narrow emission lines and take 

V0,line,l/V0j> = 0.01, k! = 1000, m = -8 (8) 
and 

foime^M.p = 100, h 2 = 1, n 2 = 0. (9) 

Shown in Fig. 1 are the \og{vf v ) — \og(v/vo tP ) curves of the fireball containing in its 
rest frame radiation the two emission lines, for various values of T. The figure shows that, 
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in the spectrum of a relativistically expanding fireball, rest frame narrow emission lines 
indeed shift to higher energy bands and are significantly broadened. Therefore, there 
would be no narrow emission lines expected in the spectrum of the object. In addition, 
a rest frame narrow emission line at high energy bands would form an up-turning tip of 
the high energy tail in the observed spectrum (an observed MeV tip can come from a 
rest frame keV line when T = 1000 or from a lOkeV line when T = 100) (see Fig. 1). 

To show the effect of the broadening, we calculate the relative width of the frequency 
range related to the FWHM of the line feature, Av FWHM /vu ne , where vu ne is the 
observed line frequency at which the peak flux of the line feature is found. Listed in 
Table 1 are the values of i>u ne as well as Avfwhm feline, coupled to various values of F, 
for the two lines. 

Table 1. List of vu ne and Avfwhm Inline for the two narrow emission lines 





Valine./ V0,p — 


0.01 




f0,(ine/f0,p 


= 100 




r 


Inline/ V0,p 


Avfwhm 


feline 


Inline 


Al/FWHM 1 Inline 


1 x 10° 


1.00 x 10" 2 


2.35 x 10" 


-2 


1.00 x 10 2 


2.35 x 10" 


2 


2 x 10° 


3.65 x 10~ 2 


1.79 x 10" 


-1 


3.65 x 10 2 


1.79 x 10" 


1 


1 x 10 1 


1.95 x 10" 1 


1.86 x 10" 


-1 


1.95 x 10 3 


1.86 x 10" 


1 


1 x 10 2 


1.96 x 10° 


1.86 x 10" 


-1 


1.96 x 10 4 


1.86 x 10" 


1 


1 x 10 3 


1.96 x 10 1 


1.86 x 10" 


-1 


1.96 x 10 5 


1.86 x 10" 


■1 


1 x 10 4 


1.96 x 10 2 


1.86 x 10" 


-1 


1.96 x 10 6 


1.86 x 10" 


■1 



We find that vu ne is a linear function of T, approximately following vu ne k, 2Tiy ,u n e, 
when T is large enough (say, T > 10). Therefore, if vu ne is identified, then one can take 
it as an indicator of the expansion speed of the fireball. It also shows that, in the case 
discussed above, where the relative width of the rest frame narrow emission lines is 0.0235 
(see the case of T = 1 in Table 1), the relative width of the line feature becomes almost 
one magnitude larger than that of the rest frame emission line when the fireball expands 
relativistically. When T is large enough, the relative width would approach an asymptotic 
value (here it is 0.186). 

As illustrated by an analytical calculation (see Appendix A) , vu ne of a very narrow 
rectangle emission, which is taken to represent a very narrow emission line, can be deter- 
mined by vune — t / o,(irie/r(l — (3). When T ^> 1, it approaches 2IV 0j im e , which is the same 
as that shown in Table 1. In this situation, the value of Av FWH M jv line can be estimated 
by 0. 162+0. 415Av 0Mne /v 0Mne . When Au 0tline /i/ 0tHne 0, Av FWHM /v Une ~ 0.162. The 
corresponding value shown in Table 1 is slightly larger than this. It might be due to the 
width of the rest frame line feature. When calculating the line described by (9) in the case 
of F = 100 by replacing = with 712 — —2, we really obtain Apfwhm feline — 0.162. 
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Fig. 2. The expected spectrum of a fireball with its rest frame radiation containing the 
two weak emission lines defined by (10) and (11), where other parameters and the symbols 
are the same as those adopted in Fig. 1. 

3. Effect on other lines 

We consider in the following the effect on other lines. 

a) Weak narrow emission lines 

We take 

vo,Une,i/vo, p = 0.01, h\ = 10, n\ = -8 (10) 
and 

vo,line,2/vo, P = 100, h 2 = 0.01, n 2 = 0. (11) 

Shown in Fig. 2 are the \og{vf v ) — \og{u/vQ iP ) curves with the two weak emission lines. 
The figure shows that, in the spectrum of the fireball considered above, rest frame weak 
narrow emission lines also shift to higher energy bands but are significantly smoothed 
and would hardly be detectable. 

b) Broad emission lines 
We take 

VQ,line,\ I 'vq, p = 0.01, hi = 10, m = -5 (12) 
and 

V0,line,2/"0,p = 100, h 2 = 0.01, n 2 = 3. (13) 

Shown in Fig. 3 are the \og(vf v ) — \og(v/v ^ p ) curves with the two broad emission 
lines. It shows that rest frame broad emission lines also shift to higher energy bands but 
their features remain almost unchanged and therefore would be easier to detect. 
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Fig. 3. The expected spectrum of a fireball with its rest frame radiation containing 
the two broad emission lines defined by (12) and (13), where other parameters and the 
symbols are the same as those adopted in Fig. 1. 

c) Narrow absorption lines 

We take 

VO,Une,l/vo,p = 0.01, h t = -1000, m = -8 (14) 
and 

^0,Kne,2/^0,p = 100, fl 2 = -1, W 2 = 0. (15) 

Shown in Fig. 4 are the \og(vf v ) — \og(v/vQ tP ) curves with the two narrow absorption 
lines. It shows that rest frame narrow absorption lines also shift to higher energy bands 
but arc significantly smoothed and would hardly be detectable. 

d) Broad absorption lines 
We take 

V0,line,l/V0,p = 0-01, h\ = -10, Tl\ = -5 (16) 
and 

V0,line,2/V0,p = 100, h 2 = -0.01, n 2 = 3. (17) 

Shown in Fig. 5 are the log(i//„) — log(v/vo lP ) curves with the two broad absorption 
lines. It shows that rest frame broad absorption lines also shift to higher energy bands 
but are significantly smoothed so that deep gaps would be significantly filled. The broad 
absorption features would remain and therefore would be detectable. 

e) Weak broad absorption lines 
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Fig. 4. The expected spectrum of a fireball with its rest frame radiation containing the 
two narrow absorption lines defined by (14) and (15), where other parameters and the 
symbols are the same as those adopted in Fig. 1. 



Fig. 5. The expected spectrum of a fireball with its rest frame radiation containing the 
two broad absorption lines defined by (16) and (17), where other parameters and the 
symbols are the same as those adopted in Fig. 1. 

We take 

^CUme,l/^0,p = 0.01, hi = -2, m = -5 (18) 
and 



vo,line,2/vo, P = 100, h 2 = -0.002, n 2 = 3. 



(19) 



Y.-P. Qin: Emission and absorption lines of gamma-ray bursts 



9 



Fig. 6. The expected spectrum of a fireball with its rest frame radiation containing the 
two weak broad absorption lines defined by (18) and (19), where other parameters and 
the symbols are the same as those adopted in Fig. 1. 

Shown in Fig. 6 are the log(^/„) — log(^/^o, p ) curves with the two weak broad ab- 
sorption lines. It shows that rest frame weak broad absorption lines also shift to higher 
energy bands but their features remain almost unchanged and would be detectable. 

f) Emission line forest 

We consider two emission line forests at lower and higher energy bands respectively. 



We take 

V0,line,l/V0,p = 0.01, hi = 1000, m = -8, (20) 

vo,lme,2/vo, P = 0.012, h 2 = 1000, n 2 = -8, (21) 

^,1^,3/^ = 0.014, /i 3 = 1000, n 3 = -8, (22) 

vo,UneA/vo, P = 0.016, hi = 1000, n 4 = -8, (23) 
and 

vo,iine,s/vo, P = 100, h 5 = 1, n 5 = 0, (24) 

vo,Unefi/vo, P = 120, he = l, n 6 = 0, (25) 

V0,line,7/v<3,p = 140, h 7 = 1, ^7 = 0, (26) 

V0,line,8/V0,p = 160, ^8 = 1, n s = 0. (27) 



Shown in Fig. 7 are the \og{vf v )— log(i//i/o, p ) curves with the two emission line forests. 
It shows that rest frame emission line forests also shift to higher energy bands but they 
are significantly smoothed and each of them forms a single broad line feature. 
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Fig. 7. The expected spectrum of a fireball with its rest frame radiation containing the 
two emission line forests defined by (20) — (27), where other parameters and the symbols 
are the same as those adopted in Fig. 1. 

g) Absorption line forest 

We consider two absorption line forests at lower and higher energy bands respectively 



and take 

v ,line,i/vo,p = 0.01, h x = -1000, n x = -8, (28) 

^(Mme,2/Vo, P = 0.012, fo 2 = -1000, n 2 = -8, (29) 

V0,line^/V0,p = 0.014, fo 3 = -1000, n 3 = -8, (30) 

vo,UneA/vo, P = 0.016, fo 4 = -1000, n 4 = -8, (31) 
and 

V , line, 5 / v , p = 100, h 5 = -1, n 5 = 0, (32) 

vo,linefi/vo,p = 120, ft 6 = -1, n 6 = 0, (33) 

V0,line,7/V0,p = 140, ft 7 = -l, n 7 = 0, (34) 

t'0,Kne,8/l / 0,p = 160, fog = -1, n 8 = 0. (35) 



Shown in Fig. 8 are the \og{vf v ) — \og{v/v 0tP ) curves with the two absorption line 
forests. It shows that rest frame absorption line forests also shift to higher energy bands 
but are significantly smoothed and would hardly be detectable. 
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Fig. 8. The expected spectrum of a fireball with its rest frame radiation containing 
the two absorption line forests defined by (28) — (35), where other parameters and the 
symbols are the same as those adopted in Fig. 1. 

4. Effect of time dependence of the line intensity or variation of the Lorentz 
factor 

Here we investigate the effect of time dependence of the line intensity and the effect of 
variation of T. 

a) The case of the line intensity varying with time 

We consider rest frame radiation containing a constant GRB form and two time 
dependent intensity narrow lines. In this situation, the formulas adopted in section 2 are 
valid when we replace hi in (7) with hi(to,e)- Suppose the line intensity is an exponential 
function of time: 

hi(t ,e) = h ifi exp(-— — — ) (to,e > *o,c), (36) 

TO 

where hifi, to, c and To are constants. As shown in Paper I, t .e and t can be related by 

_ t-t c -D/c+(R c /c) cos 6 

to ' e r(i-/?cos0) + io ' c ' (37) 

where t c and R c are constants. We take 

r = 100R c /c. (38) 
Therefore, 

h( , , . r C(t l: ) ' g + cosg 1 r t-t c - D/c+(R c /c) COS g 

h tM = ^oexp[- ioor(i /3cosg) ] [ r{1 _ poQBe) > 0]. (39) 

We consider different observation times and take 

ti = —+t c , (40) 
c 
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Fig. 9. The expected spectrum of a fireball containing in its rest frame radiation the two 
time dependent intensity emission lines [confined by (39)] for observation time t\ [defined 
by (40)], where we take 2ttIqvq^R 2 (t\) / D 2 = 1. The solid lines from the bottom to the 
top correspond to Y = 1,2, 10, 100, 1000, 10000, respectively, and the dotted lines (some 
are totally overlapped by the corresponding solid lines) are those presented in Fig. 1. 

t2 = ^±R +tc (41) 
c 

and 

WR C + D 

h = + t c . (42) 

c 

Once more we adopt (8) and (9), where hi stands for hi t o (i = 1,2). 

The plots of \og(i>f v ) — log(i//i/o iP ) for the fireball containing in its rest frame radiation 
the two time dependent intensity emission lines for the three observation times are shown 
by Figs. 9-11 respectively, where various values of T are considered for each figure. For 
the sake of comparison, we take 2irI v 0p R 2 (ti) / D 2 — 1, 2irlQVQ, p R 2 (t2) / D 2 — 1 and 
2-kIqVq^R 2 (t%) I D 2 = 1 for the three figures, and the curves shown in Fig. 1. are also 
plotted. 

Similar to that shown in Fig. 1, in the spectrum of a relativistically expanding fireball, 
rest frame dimming narrow emission lines would also shift to higher energy bands and 
would be significantly broadened. However, for fireballs with a very large Lorentz factor 
(say, T > 1000), the line feature would disappear. As time progresses, the feature would 
eventually disappear and this would happen at a much earlier time for those fireballs 
with a relatively larger value of the Lorentz factor. 

b) The case of variation of the Lorentz factor 

Here we study the effect of variation of the Lorentz factor on the spectrum of fireballs. 
The expected spectrum of fireballs influenced by this variation is studied in Appendix 
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Fig. 10. The expected spectrum of a fireball containing in its rest frame radiation the 
two time dependent intensity emission lines for observation time f 2 [defined by (41)], 
where we take 2-kIqvq^R 2 fa) / D 2 = 1. Other parameters and the symbols are the same 
as those adopted in Fig. 9. 



Fig. 11. The expected spectrum of a fireball containing in its rest frame radiation the 
two time dependent intensity emission lines for observation time t 3 [defined by (42)], 
where we take 2-kIqv ^ p R 2 fa) / D 2 = 1. Other parameters and the symbols are the same 
as those adopted in Fig. 9. 

B. The effect would be shown more clearly when the influenced spectrum is compared 
with the previously discussed one. Therefore, we consider the same rest frame constant 
radiation discussed in section 2, with the variation of the Lorentz factor being the only 
difference. To employ what we derived in Appendix B, we assume that the radiation is 
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independent of direction, Io,v{to.e,^o.e,&) = Io,v (to, 0,^0,8), and there is no limit for the 
radiated frequency and the radiation concerned lasts a sufficiently large interval of time 
so that the integral limits of (B.15) are 9 mm — and # max = n/2 (see Paper I for a 
detailed discussion). Applying (3), (B.7) and (B.16), we obtain from (B.15) that 

vm = —Jo fWeW^m^w cos d ' (43) 

where we take Io(to,e) — lo, and go, 1/(^0,0) is determined by (6). Frequencies vo,g and v 
are related by the Doppler effect [see (B.7)]. 

It is known that, in its deceleration phase, the Lorentz factor of a fireball would 
develop with time following i~ 3 / 8 in the case of the adiabatic hydrodynamics while fol- 
lowing i~ 3 / 7 in the case of the radiative hydrodynamics (Piran 1999). To illustrate the 
effect, we consider the latter case and assume 

T(tg) = m&x{ktg 3/7 , 1} (tg > t c ), (44) 

where k is a constant. Here we assign t c to be the time when the deceleration starts and 
assign T = 1 when the deceleration as well as the expansion stop. Suppose T c > 1 is the 
value of the Lorentz factor of the fireball at time t c , then one obtains 

k = T c tl' 7 . (45) 

When T{te) is known, the relation between the coordinate time, tg, of the differential 
surface of the fireball concerned and the observation time, t, will be well established, and 
the radius of the fireball as a function of time will be determined (see Appendix B). From 
(B.9) we get 

U [ 1^ \F- ~-, jJ-Vr -^ dte] ™se = t+^cos6-°. (46) 
Jt c y (max{fci e ' ,l\Y c c 

We find that, once t is assigned, tg will be determined by (46). Then T(tg) will be 
determined by (44). Once both t and tg are known, R(tg) will be easily determined. We 
get from (B.8) and (46) that 



= ° ° {t te) =R c + cf\ ll- 1 dtg. (47) 

cos 9 J tc y (max{fci~ 3/7 , l}) 2 

To show the effect, different observation times should be concerned. Let 

tic 

Then 

D+pR c 

t = t c + — . (49) 

c 

We then come to 




tg-[l Jl \r 7 dt ]cos8 = t c + — ( P + cos9) (50) 

' 1 (max{fct; 3/7 ,l}) 2 c 
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Fig. 12. The expected spectrum of a fireball with a variable Lorentz factor, where the 
initial Lorentz factor is r c = 2 and the radiation concerned is that discussed in section 
2. We take R c /c = 1, 2irI v a , p R 2 c /D 2 = 1, and k = 10k (p = 10, T c = 2,R c /c = 1). 
The solid lines from the bottom to the top correspond to p = 0, 10, respectively, and 
the dotted lines (some are totally overlapped by the corresponding solid lines) are those 
presented in Fig. 1. 

and 

To focus on the effect of variation of the Lorentz factor, we require that, all photons 
observed at the assigned time must be those emitted after the deceleration starts and 
before the expansion stops. The first requirement leads to p > (see Appendix C). To 
meet the second requirement, k cannot be arbitrarily taken as long as p, T c and R c are 
fixed. For given values of p, T c and R c , there will a lower limit of k (see Appendix C), 
i.e., 

k > k , (52) 
with 

fco = r c [ ]3A (53) 

c (i-/3 c )(rI /3 -i) 

where C = ^/Tf^l/T c . 

To show the development of the spectrum, we consider two different observation 
times and take p = (0, 10). For the sake of comparison, we adopt R c /c = 1 and take 
2irI v ^ p R 2 / D 2 = 1 and k = !Qk {p = 10,r c ,i? c ) for all situations concerned. 

The results are displayed in Figs. 12-16. All these figures show that the flux observed 
at a relatively later time would be larger than that observed at an earlier time, suggesting 
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Fig. 13. The expected spectrum of a fireball with a variable Lorentz factor, where L c = 10 
and k = I0ko(p = 10, T c = 10, R c /c = 1). Other parameters and the symbols are the 
same as those adopted in Fig. 12. 



Fig. 14. The expected spectrum of a fireball with a variable Lorentz factor, where T c = 
100 and k — 10k (p = 10, L c = 100, R c /c = 1). Other parameters and the symbols are 
the same as those adopted in Fig. 12. 



that, in the situation studied here, the enhancement of the radius of the fireball plays an 
important role in the magnitude of the flux. 

Shown in Fig. 12 we find that, when T c is very small, as the variation of the Lorentz 
factor is then not obvious, the effect of deceleration would be insignificant. The positions 
of emission lines would remain almost unchanged. 
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Fig. 15. The expected spectrum of a fireball with a variable Lorentz factor, where T c = 
1000 and k = 10fco(p = 10, r c = 1000, R c /c — 1). Other parameters and the symbols are 
the same as those adopted in Fig. 12. 



Fig. 16. The expected spectrum of a fireball with a variable Lorentz factor, where T c = 
10000 and k = I0k a (p = 10, T c = 10000, R c /c = 1). Other parameters and the symbols 
are the same as those adopted in Fig. 12. 



Figs. 14-16 reveal that the flux affected by the deceleration would become much 
smaller than the constant radiation flux when the initial Lorentz factor is large enough 
(say, r c > 100), and the change of the positions of emission lines due to the decrease 
of the Lorentz factor would be quite significant (the positions of the lines would shift to 
lower energy bands as time progresses). 
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Fig. 17. The expected spectrum of a fireball with a variable Lorentz factor, where we 
take p = 10 and k = 10k n (p — 10, T Cl R c /c = 1) and consider different values of the 
initial Lorentz factor L c . The solid lines from the bottom to the top correspond to T c = 
10, 20, 50, 100, 200, 500, 1000, respectively, and the dotted lines from the bottom to the 
top are the T = 10 and V = 100 lines presented in Fig. 1. Other parameters are the same 
as those adopted in Fig. 12. 



Fig. 18. The expected spectrum of a fireball with a variable Lorentz factor, where we take 
r c = 300 and k = 10fco(j> = 10, r c = 300, R c /c — 1) and consider different observation 
times, which are represented by different values of p, respectively. The solid lines from 
the bottom to the top correspond to p = 0, 1, 3, 10, 100, 1000, 2000, respectively, and the 
dotted lines from the bottom to the top are the L = 100 and L = 1000 lines presented in 
Fig. 1. Other parameters are the same as those adopted in Fig. 12. 
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We find in Fig. 14 an unfamiliar form of emission lines. To view this phenomenon in 
more detail, see Fig. 17, where more lines with different values of T c are presented (a 
detailed study of this issue will be made later). It shows that, for the assigned observation 
time, the emission lines concerned would be sharp on both edges in some cases where the 
initial Lorcntz factor is carefully chosen. To find out if this phenomenon depends also on 
the observation time, we plot Fig. 18, where different observation times are considered 
for a single initial Lorentz factor. One finds that the shape of the lines is indeed changed 
with time. During some period, it would be sharp on both edges. Meanwhile, the flux 
keeps being enhanced (at least during the epoch of the deceleration phase concerned). 
This would become an interesting expected observational characteristic. 

5. Discussion and Conclusions 

In this paper, we investigate how the Doppler effect in the fireball framework affects the 
observed spectrum of the object with various rest frame emission or absorption lines. 

In section 2, we study the effect on narrow emission lines. The study reveals that, 
influenced by the Doppler effect in the fireball framework, rest frame narrow emission lines 
would significantly shift to higher energy bands and would be significantly broadened. 
A rest frame narrow emission line at high energy bands could form an up-turning tip of 
the high energy tail in the observed spectrum. It shows that, when T is large enough, the 
observed line frequency nine and the rest frame line frequency fo,iine would be related by 
vn ne « 2IVo,;me- It also reveals that the relative width of the line feature would approach 
an asymptotic value when T 3> 1. In the case where the relative width of the rest frame 
narrow emission lines is 0.0235, the observed Avfwhm /^Une would become almost one 
magnitude larger than that of the rest frame emission line. When the rest frame line is 
narrow enough, the relative width would approach 0.162. 

In section 3, the effects on other lines are investigated. The study shows that in the 
spectrum of a relativistically expanding fireball, any rest frame lines would shift to higher 
energy bands and would be significantly smoothed. Rest frame weak narrow emission lines 
as well as narrow absorption lines and absorption line forests would be smoothed and 
would hardly be detectable. Meanwhile, the features of rest frame broad emission lines 
as well as both strong and weak broad absorption lines would remain almost unchanged 
and therefore would be easier to detect. Deep gaps caused by rest frame broad absorption 
lines would be significantly filled. In addition, each of the rest frame emission line forests 
would form a single broad line feature. 

In section 4, we investigate the effect of time dependence of the line intensity and 
the effect of variation of the Lorentz factor. The study reveals that, in the spectrum 
of a relativistically expanding fireball, rest frame dimming narrow emission lines would 
also shift to higher energy bands and would be significantly broadened. However, when 
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the Lorentz factor is very large (say, T > 1000), the line feature would disappear. As 
time progresses, the feature would eventually disappear and this would occur at a much 
earlier time for those fireballs with a larger value of T. Considering the variation of 
the Lorentz factor of fireballs we find that the flux observed at a relatively later time 
would be larger than that observed at an earlier time, suggesting that, in the situation 
studied here, the enhancement of the radius of the fireball plays an important role in 
determining the magnitude of the flux. When the initial Lorentz factor is very small, the 
effect of deceleration would be insignificant, and the positions of emission lines would 
remain almost unchanged. When the initial Lorentz factor is large enough, due to the 
deceleration, the flux would become much smaller than the constant radiation flux and 
the change of the positions of emission lines would be obvious. As time progresses, the 
positions would shift to lower energy bands. In addition, we find an unfamiliar form of 
emission lines that is sharp on both edges. Further study reveals that this phenomenon 
depends not only on the initial Lorentz factor but also on the observation time. For 
some initial Lorentz factors, the shape of lines changes with time. This might become an 
interesting observational characteristic expected in GRBs. 

One can conclude that, for a rclativistically expanding fireball, there would be no 
narrow emission lines expected in its spectrum, and rest frame narrow absorption lines 
as well as rest frame emission or absorption line forests would not be detectable. Since 
vu ne is a linear function of T, features of lines could serve as an indicator of the expansion 
speed of fireballs as long as they are identified. In the phase of deceleration of the Lorentz 
factor, one can expect a change of the shape of emission lines of fireballs with time. 

If there are some high energy emission lines such as the 6AkeV line and the bllkeV 
annihilation line in the outer shell of fireballs, we could expect them to appear at around 
O.lMeV and lOMeV respectively when T ~ 10, or lMeV and lOOMeV when T - 100, 
or lOMeV and lOOOMe^ when T ~ 1000. Regarding the previous report of the detection 
of emission features near 400keV (Mazcts et al. 1980), we would prefer to interpret them 
as the blue-shifted QAkeV line corresponding to the fireball expanding with T ~ 30. Due 
to the great blue shift effect, the conventional redshift interpretation of the line seems 
unlikely. We propose that, in detecting line features, the blue shift and broadening effects 
due to the expansion of the fireball should be taken into account. 

In addition to the conventional fireball model, the method can also be applied to the 
cannonball model (Dado et al. 2002a, 2002b). To explain the X-ray line features observed 
in some GRBs as being Doppler-boosted, Dado et al. (2002c) employed a cannonball 
model for which the radius of the cannonball was treated as a constant and the direction 
of the motion of the cannonball was assumed to be unchanged. In this situation, some 
formulas in Paper I can be simply applied. Taking i ,e — *o, c as the proper emission time 
to and t — t c — D/c+ (R c /c) cos9 as the observation time t we would get from (A. 8) of 
Paper I that t — S(t)t, where S(t) = 1/T(t)(l — f3(t)cos9) is the Doppler factor of the 
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cannonball at time t. The observed area of the cannonball is S(9,T)R 2 , where S(9,T) 
is a function of the Lorentz factor, which would be it when the contraction effect due 
to the relative motion is ignored. Replacing R 2 (ts) cos 9 sin 9d9d(p with S(9,T)R 2 we get 
from (A. 11) of Paper I that dE = /S.s bdvdtS{9,Y)R 2 I u {v,t) / D 2 , where the cosmological 
effect is ignored. Absorbing S{9, V) into the normalized coefficient we get F b s (i/,t) = 
R 2 I v (v, t)/D 2 . Applying the Dopplcr effect and replacing I v {v, t) with <5 3 (t)/o,!/(^o, to) = 
S 3 (t)FcB[i / /5(t), 6(t)t] we come to equation (3) of Dado et al. (2002c) and the Dopplcr 
effect itself leads to equation (6) there (when the cosmological effect is ignored); based 
on this the spectral data (including the X-ray line features) of some GRBs were fitted 
(Dado ct al. 2002c). 

Acknowlegements. It is my great pleasure to thank Prof. M. J. Rees for his helpful 
suggestions and discussion. This work was supported by the Special Funds for Major State 
Basic Research Projects ("973") and National Natural Science Foundation of China (No. 
10273019). 

Appendix A: Analytical study of AfFWHAi/^Une f° r ver Y narrow emission 
lines 

Here we calculate IS.ufwhm Inline analytically for very narrow emission lines when r>l. 

The flux expected from a fireball with a definite value of the Lorentz factor is (see 
Paper I) 



where R{t) is determined by (2); t, v, D, 9, vo.e, io,e, and Io^(to,e, vofi) are explained in 
section 2. As discussed in section 2, we consider a constant radiation and get the following 
by applying (3): 



To make an analytical study of Avfwhm jviine for very narrow emission lines, one 
would prefer a simple form of emission, which can be calculated analytically and can 
represent an emission line in extreme cases. 

Let us consider a very narrow rectangle of emission. The rest frame radiation is 
assumed to be 

9o,v{va,e) = go, v ,c(vo,e) + e 0tV {vo#), (A. 3) 

where, while go,v,c{ v o,e) describes the relative intensity of a continuum (see section 2), 
eo, i/(^o, e) stands for the relative intensity of an emission line which is 

eo,i/(fO,0) = h inline < Vq,8 < VQ.line + Al/Q.Mne), (A. 4) 




(A.I) 




(A.2) 
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where ft. is a constant. 

To calculate Aufwhm /vime, one needs only to deal with the radiation of the line: 

fv,line{t) = £)2^3 " &v ( A -5) 

with 

nax cos 6» sin .9 , k ^ 

f3cos9) 5 y ' 

According to the Doppler effect, for certain values of vo,e and T, the shifted frequency 

v is confined in the range VQ.ejY < v < ^o,e/r(l — /?) for an expanding fireball. Thus, 

the valid range of frequency in (A. 6) is 

line / \ rr\ 

— r(i-/3) • (AJ) 

The integral limits # m i n and 9 max in (A. 6) can be determined by (Paper I) 

flmi„ = cos- 1 (min{l, i(l - ^ip)}) 

(0<i(l-^);0</3) 

and 

C = cos" 1 (max{0, 1(1 - )}) 

(i (1 _te+£te)<i ;0</ 3). 
Integrating (A. 6) we get 



(A.E 



(A.9) 



_ A/lr i 1 l 

/5 2l 4 L (l-/3min{l i (1- )})4 (1-/3 max { , "".»"^ + ^"" dHl£)})4 ^ 

■ if i i n (A 10) 

3[ (l-/?max{0,4(l- t '"^— tr - 2 ^)}) 3 (l-/3min{l i(l-lH^)})3^ 



r — " — r(i-/3) 



v r - " - rn-m J- 



Let us calculate de v jdv in different frequency ranges: range I, ( ""'{!" < f < 

.0,H„e+A,o,H„ e ); range ( -O.H„e+A, ,H„ e < y < _^ ); ^ range ( < y < 

ly °''"r("i-^)'''" e )■ Here, we assume 

AU0 ' lme « 1 and l«r, (A.ll) 

which leads to 

^- +A ^-«T^- (A.12) 
Applying (A.12), we find de v /dv > in both ranges I and II, and de v jdv < in range 
III. This leads to 

4 

OH "0,lme 1 3 

_ h f [ (^O.iirve+At/Q.H^^J , V Q .line -,1 /A 1 Q\ 

e,, max - 3/?2(1 _ ^3 i 4(1 _ /?) + ^ ^ + Al ^ iMne)3 A >. l A '^ 

which is the maximum value of e„ in the whole frequency range. 

The frequency at which e„ jmax is found is identified as the observed line frequency 
vu ne which is determined by 
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When r 3> 1, it would approach uu ne ~ 2FVo.Zme- 

With (A. 13), we first calculate vfwhm in region II, vfwhmji, and obtain 

„4 4 4 4 

1 V FWHM,II f-i u O,lin<: 



ivl Une L " (^0,H„e + A^,H„e) 4j 8(1-/3)4 t 1 , H„ e + A^O , H„e ) 4 J 

_ r ^FiyffM.j-J r-i ^O.linc 1 _ 1 r-i ^OMne 1 

_ 3< H „ e 1 (^l,»e+A"0 |1 ,„e) !| 6(1-/3)3 L 1 („„ , , ~ " + A U 0> J~ ) 3 J ' 

The solution of (A. 15) can be approximated by 

Tvfwhmji „, [l + W 1/12 -l)] 1/3 u | 2(2V 12 -1) A^ m 

HUine " 21/3(1-/3) l + 3(2 1 / 12 — 1) fo,Hne ^ ' 

Second, we calculate vfwhm in region III, vfwhmjii, and get 



(A.15) 



4 4 1 + (1+ T°'""° ) 4 

' FW H M ,1 1 1 v ,y O,HTie _ 



Kiine 8(l-/3) 4 



(A.17) 



The solution of (A.17) can be approximated by 

TVFWHMJII 1 /r.l/3 , n-2/3 ^O.Une ^ 



1 Al/QJ, 



(2 1 / 3 + 2- 2 / 3 °'' me ) = 2 (A.18) 



We therefore obtain 

~ 1 _ r l+3(2 1 / 12 -l) 1 i/ 3 , r 1 _ 2 2 / 3 (2 1 / 12 -l) , 
v line — 1 L 2 J T l2 [l+3(2i/i2-l)]2/3 J ^ 



~ 0.162 + 0.415^-' 



(A.19) 



fO.line 

When Avo,line/vo,line — > 0, One will get Al/FWHAl/viine — 0.162. 

Appendix B: Expected flux of a fireball with a variable Lorentz factor 

We examine a fireball expanding at a variable velocity v = (3c and employ the same 
coordinate system adopted in section 2. Consider radiation from a rest frame differential 
surface, c?so,0.</" of the fireball at proper time i ,e- Let ds e ^ be the corresponding differ- 
ential surface resting on the observer's framework, coinciding with dso,e, v at to^g, and let 
tg be the corresponding coordinate time of to.e- Obviously, tg = tg(tofi). 
According to the theory of special relativity, tg and to,e are related by 

dt e = T(t e )dt ,8 = r (to,e)dt ,9, (B.l) 

where T(tg) is the Lorentz factor of the fireball at tg and r (io,#) = ^\pe{to,e)\- Integrating 
(B.l) yields 

tg = I T (t ,e)dt ,e + t c or t ,g = / T7TT dt( > + *o,c (B.2) 

Jto.c Jt c r W 

here we assign tg = t c when t ,6" = to,c- 
The area of dsg jlp is 

dse iV = i? 2 (t e ) sm0d0d<p, (B.3) 

where i?(ig) is the radius of the fireball at tg. Suppose the radius develops as 

dR(tg) = Cp(tg)dtg. (B.4) 
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Integrating (B.4) leads to 

R(t e ) = c J " /3(tg)dt g + R c , (B.5) 

where R c is the radius at time tg = t c . As assigned above, t$ and to t $ correspond to the 
same moment. Therefore the radius can be expressed as 

rta(to,e) r^o,e 

Ro(t Oi0 ) = R[t (t Oi0 )] = c / 0(to)dt e + Rc = c / 0o(to,e)r o (t o ,g)dt Oi g + R c , (B.6) 

Jt c Jt ,c 

where (B.l) is applied, and /3b(*o,e) = P[te(to,e)]- 

Let us consider an observation within small intervals t — t+dt and v — v+dv carried 
out by an observer with a detector As b at a distance D (D is the distance between the 
observer and the center of the fireball), where D 3> R(tg). Suppose radiation from dso.g yip 
reaching the observer within the above observation intervals is emitted within the proper 
time interval io,0 — to,e + dtoj and the rest frame frequency interval 1/0,0 — ^0,0 + dv .g. 
According to the Doppler effect, v and 1/0,0 are related by 

^0,0 ^0,0 



" r(i o )[l-/3(i 9 )cos0] r (t , 9 )[l-A)(to,o)cos^]' 

Considering the travelling of light from the fireball to the observer we get 



(B.7) 



c(t-tg) = D - R(tg) cos 6 (B.8) 
(here the cosmological effect is ignored). Combining (B.5) and (B.8) we get 

I3{tg)dtg] COS 6 — tg = D / c — t — (R c /c) COS 6, (B.9) 

while combining (B.2), (B.6) and (B.8) we obtain 

rto,e /■*<), e 

[/ 0o(to,e)r o (t o ,e)dto,e]cos6-[ r (t ,e)dt 0i g] = D/c-t- (R c /c) cos6> + i c .(B.10) 

Once (3(tg) or /3o(to,e) is known, tg and t o ,0 as functions of 9 and t can be determined 
by (B.9) and (B.10), respectively. Thus the radius of the fireball can be determined by 
(B.5) or (B.6). 

Suppose photons, which are emitted from dso,e yip within the proper time interval to,e 
— t ,g + dt ,g and then reach the observer within t — t + dt, pass through dsg_ v within 
the coordinate time interval tg — tg + dtg :S . Since both dsg ylfi and the observer rest in the 
same framework, dtg, s = dt. Of course, the frequency interval for the photons measured 
by both dsg tV and the observer must be the same, which is taken as v — v +du. In views 
of dsg, v , the amount of energy emitted from dso : e yip within to,e — to,0 + dto,e and ^0,0 - 
1/0,0 + dvQ.g towards the observer would be 

dEg tV = I lJ (tg 1 v,9)cos9dsg yV dujdvdt, (B-ll) 

where I v (tg,v, 9) is the intensity of radiation measured by dsg :V , which is allowed to be a 
function of direction, and du is the solid angle of As b with respect to the fireball, which 
is 

*» = ^- (B.12) 
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Thus, 

, v, 9) cos 9 sin 9d9dip 

d.be, v = jy^ , (B.13) 

where (B.3) and (B.12) are applied. 

The total amount of energy emitted from the whole fireball surface detected by the 
observer within the above observation intervals is an integral of dEg tV over that area, 
which is 

dE = 2nAs ^dt [ 6 ™ x R 2( tMtf) 0) cos sin 0d0, (B.14) 

^ $min 

where 9 m i n and 9 max are determined by the fireball surface itself together with the emitted 
ranges of i ,e and is ,e- Thus, the expected flux would be 
2tt /" emax 

f^ = ly2 R 2 (t e )I»(tg,is,9)cos9sm9d0. (B.15) 

J ^min 

It is well known that the observer frame intensity I v {te,v,6) is related to the rest 
frame intensity Io. v (to.e,i'o,e,9) by 

I„(t e ,is,9) = (—fl Q ^t .o,v , e ,9). (B.16) 
The flux then can be written as 

t M _ 2?r / >8max Rl{t ,e)Io^{t , e ,vo.e,9) cosflsin6> Jfl 

Mj ~W 9min r§(t 0l9 )[i-A)(to, e )cos(?]3 ^ i7j 

where (B.6) and (B.7) are applied. 

The range of 9 of the visible fireball surface is 

< 9 < tt/2. (B.18) 

Within this range, suppose the emitted ranges of t ,e and ^ ,e constrain by 

#t,min < 9 < 9 t , maK (B.19) 

and 

#;/,min < 9 < ^^max, (B.20) 

respectively. Then when the following condition 

max{6> timin , 9^ min } < min{6> t , max , 9 Utmax _} (B.21) 
is satisfied, ^min and 9 max would be obtained by 

#min = niax{0 t nl i n , ^min} (B.22) 

and 

u max 

= min{# 4 
respectively. 
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Let the emitted ranges of to,e and vo,e be 

io.min < to,e < ^0,max (B.24) 

and 

^(^min .max? 

(B.25) 

respectively. For > 0, one can obtain the following from (B.7) and (B.25): 
^,min(io,e) = cos ^ (min{l, ±(1 - ^r-)}^ 



(0<i(l-^);0</3 
and 

Ov, max (to,e) = cos" 1 fmax{0, - ] ^f L )} 



(B.26) 



(B.27) 

(i(i-^)<l;0</3), 

where = 0{t e ) = o {t o ,e) and T = T(t g ) = T (t ,e)- From (B.10) and (B.24) one can 
also determine #t im m and #t,max once 0{t$) or 0o(to,e) is known. 

Appendix C: Lower limit of the coefficient of deceleration 

We consider in section 4 the case of the radiative hydrodynamics for which, when decel- 
erated, the development of the Lorentz factor of a fireball follows that shown by (44). 
The coordinate time, tg, of the differential surface of the fireball concerned and the ob- 
servation time, t, are related by (B.9). Since 0(t$) < C , where C = y/T^ — 1/T C , we 
obtain from (B.9) that 

t - R + (Ec _ t p ) COS 

to < — c —r^ — Y 1 — • C.l 

1 — 0c cos 

Applying (49) we get 

R c p + cos 9 ,_, . 

te<t c + ^ f— -. C.2 

c 1 — 0c COS 

To focus on the effect of variation of the Lorentz factor, we require that all photons 
observed at the assigned time must be those emitted after the deceleration starts and 
before the expansion stops, i.e., 

t e > t c (C.3) 
and 

T(t e ) > 1. (C.4) 

Combining (C.2) and (C.3) we get 

p>0. (C.5) 

One can verify by applying (44) and (45) that, when 

k^>t c + ^-^pl-, (C.6) 
c 1 — p c cos 
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then (C.4) would be satisfied. Taking 9 = we get 

k r /3>tc+ R 1 1+P (C7) 

C 1 jj c 

Applying (45) we arrive at 

fc>r c [ R ^ 1+ 7 f 3 (C.8) 
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